Spin and interaction effects in quantum dots: a Hartree-Fock-Koopmans approach 
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We use a Hartree-Fock-Koopmans approach to study spin 
and interaction effects in a diffusive or chaotic quantum dot. 
In particular, we derive the statistics of the spacings between 
successive Coulomb-blockade peaks. We include fluctuations 
of the matrix elements of the two-body screened interaction, 
surface-charge potential, and confining potential to leading 
order in the inverse Thouless conductance. The calculated 
peak-spacing distribution is compared with experimental re- 
sults. 

The traditional description of Coulomb blockade in 
quantum dots has been the constant-interaction (CI) 
model, in which the electrons occupy single-particle lev- 
els in a confining potential, and the interaction is taken 
as a constant charging energy. In dots with chaotic dy- 
namics, the fluctuations of the single-particle levels and 
wave functions can be described by random-matrix the- 
ory (RMT). The CI plus RMT model was successful in 
describing some of the observed statistical properties of 
such dots, e.g., the conductance peak-height distribution 
[0,1). However, several experiments have demonstrated 
that other statistics, such as the distribution of the spac- 
ing between Coulomb-blockade peaks, are affected by 
electron-electron interactions 

At low temperatures the peak spacing is given by the 
second-order difference of the ground-state energy versus 
the number of electrons. In the CI model and for spin- 
degenerate levels, the peak-spacing distribution is bi- 
modal, i.e., a superposition of a 8 function and a (shifted) 
Wigner-Dyson distribution. However, none of the mea- 
sured distributions are bimodal and they all deviate from 
Wigner-Dyson statistics J|,||,|]J|. 

In the spinless case, exact diagonalization for a small 
number of electrons |j, Hartree-Fock (HF) and den- 
sity functional Q calculations, as well as a random in- 
teraction matrix model || explained the deviation from 
Wigner-Dyson statistics as an interaction effect. Spin 
degrees of freedom were included using exact diagonal- 
ization for dots with a small number of electrons and 
small values of the Thouless conductance g 
effects were studied in the limit g — > oo [p~T|. 
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Spin 
us- 
ing the so-called universal Hamiltonian [[I4| , [l5| . Inter- 
action effects (in the presence of a magnetic field) were 
included using a Strutinsky approach for quantum dots 
p2| . The resulting distributions showed qualitative dif- 
ferences when compared with experiments. Temperature 
effects were shown to be important at temperatures as 
low as T ~ 0.1 — 0.2A |ll| (A is the mean spacing be- 
tween spin-degenerate levels) . The spacing statistics were 



also studied in a spin density functional theory for dots 
with ~ 10 electrons 

Here we develop a theory that includes spin and inter- 
action effects and is based on a HF-Koopmans approach. 
The theory is generic and does not require the actual 
solution of the HF equations. Rather than using the 
non-interacting basis, we choose as a reference state the 
n-electron dot with spin S = (n even) and work in its 
HF basis. We then consider the addition energies and the 
energy differences between various spin configurations in 
Koopmans' limit where the single-particle wave functions 
do not change ]r^| . Previously, Koopmans' approach was 
discussed for dots with spinless electrons [Q and was used 
to study spectral scrambling IpTsfl . Here we derive the 
generic statistics of the spin and peak spacings, assuming 
that the HF levels of the reference state satisfy statistics 
typical of a diffusive or chaotic dot. The theory is valid 
for large g, and for g — ► oo it reduces to the universal 
Hamiltonian |14|] . For finite g, we include fluctuations 
of the diagonal matrix elements to leading order in 1/g. 
Compared with the approach of Ref. jl2| , our theory re- 
quires the statistics of only a few levels around the Fermi 
energy, and some of its results are qualitatively different. 

The Hamiltonian of the quantum dot in the "disorder" 
basis \i a) — a\ a \0) (i denotes a spatial orbit and a = ±1 
is the spin variable) is given by 



where ej are the single-particle energies and Uij-ki are 
the (spin-independent) matrix elements of the Coulomb 
interaction. The Hamiltonian ([!]) can be solved in the HF 
approximation. For each value of the spin projection S z , 
we use Slater determinants with n + (n_ ) spin up (down) 

(n) 

orbitals. Usually, the HF single-particle energies e a J and 
orbitals cj) a(T depend on the spin a. However for even n, 
the HF equations have a solution where e« and (f> a are 
independent of a, and the lowest n/2 levels are doubly 
occupied. Such a Slater determinant has good S = 0. 
We choose this solution as our reference state, and work 
in its HF basis \a). This S — state is an eigenstate of 
the following diagonal many-particle Hamiltonian 
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where n a<J is the number operator of the state aa, and 
v a/ 3 = v al 3- a i3, u^g = v a /3 ; p a , and w^g = u q/3 - are di- 
agonal, exchange and antisymmetrized matrix elements, 
respectively. 

The Hamiltonian (||) also has eigenstates with 5 ^= 
0. We label by a — the last occupied level of the 
5 = state. The lowest energy state for each spin S is 
obtained by promoting S spin down electrons from a = 
0, . . . , — (S — 1) to spin up electrons in a — 1, . . . , S. The 
resulting Slater determinant describes a maximal spin 
projection state S z = S and has good spin S. Using 
(g), the energy difference 6E n (S) = E n (S) - E n (S = 0) 
between the lowest states with spin S and spin S = can 

(n) 

be written in terms of e a and a few matrix elements. For 
example 



SE n (S = 1) = (4 



(n) » 
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The spin S n of the ground state of the n-electron dot is 
determined by minimizing 6E n (S). 

Similarly, the ground-state spin S n+ i of the dot with 
n + 1 electrons can be determined from the energy dif- 
ferences 6E n+1 (S) = E n+1 (S) - E n+1 (S = 1/2) for half- 
integer S. Assuming Koopmans' limit, we can express 
5E n+ i(S) in terms of the HF levels and matrix elements 
of the reference state (n, S = 0). For example 



SE n+1 (S = 3/2) = (4 n) - e<j n) ) - v w - v 20 + 



(4) 



The addition energy /x n +i = E gs (n + 1) — E gs (n) is 

Hn+i = Hn+i{0 -> 1/2) + 5-B„+i(5„ + i) - SE n (S n ) , (5) 

where in Koopmans' limit /i n+ i(0 — ► 1/2) = e^. 

The spacing A2 between successive peaks is given by 
the difference in addition energies. We have to distin- 
guish between even-odd-even ("odd") and odd-even-odd 
("even") transitions (in particle number). We consider 
here the odd transition n — >n + l^n + 2(n even), 
for which A2 = fJ. n +2 — Mn+i (similar results can be 
derived for the even transition [pj[). (J, n +i is given by 
(|J) and /i n +2 is calculated from /i n +2 = /i„+2(l/2 — > 
0) + ^„ +2 (5„ +2 ) - <5^„ + i(5„+i) with M n +2(l/2 -» 0) = 
^ (n) ^ <5^„+i(5) is given by, e.g., Eq. (§, while 



1>H. 



5-B n+ 2(S') is calculated from, e.g. 



SE n+2 {S=l) = (e 



( n )\ 1 
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To describe the statistics of A2 it is necessary to model 
the fluctuations of the HF levels and wave functions of 
the (n, S = 0) dot. The spectrum is assumed to follow 
RMT within g levels around the Fermi energy, and the 
matrix elements are uncorrelated from the single-particle 
spectrum. An exception is the gap — e o ■ We find 
its statistics by comparing the single-particle spectrum 
of the (n + 2)-electron dot with the spectrum of the n- 
electron dot (both in their S = configuration). We have 



the relation e\ 
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which the levels £q™ +2 ' ) and e^ n+2 - ) are both doubly occu- 
pied, and thus their spacing should follow Wigner-Dyson 
statistics. The gap distribution is then a convolution of 
a Wigner-Dyson distribution with a Gaussian describing 
the distribution of Vq X + V01 — wn (see below). 

We apply our HF-Koopmans approach in a restricted 
single-particle space of ~ g levels around the Fermi en- 
ergy. The long-range bare Coulomb interaction should 
then be replaced by an effective interaction. In the limit 
r s « 1 wc employ an effective RPA interaction cal- 
culated by excluding particlc-holc transitions within the 
above strip of ~ g levels [^5). This effective interaction 
is v(ri,r 2 ) = e 2 /C + v K (r ll r 2 ) + V(ri) + V(r 2 ), where 
v K (ri,r 2 ) is a two-body screened interaction in the dot, 
and V(r) is a one-body potential generated by the accu- 
mulation of surface charge in the finite dot [^FJ . 

We decompose the interaction in (^) into its average 
and fluctuating parts, and denote by Ud = v a p and 
J s — v^a {a ^ (3) the average values of the direct and ex- 
change interaction, respectively. The average interaction 
is invariant under a change of the single-particle basis 
when v aa = Ud + J s and can be written in terms of the 
number operator h and total spin S [Ell 
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The average interaction is determined by the spatial 
correlations of the single-particle wave functions p0|,22|. 
To leading order in 1/g, J s = (A/2)(l + 2f3~ 1 b 1 /g) [g3|, 
where 61 is a geometry-dependent coefficient determined 
from b\jg = J drll(r, r)/A. Here n(ri,r2) is the diffu- 
son propagator and A is the area of the dot. 
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FIG. 1. Statistics in the absence of fluctuations of the in- 
teraction matrix elements. The probabilities of various spin 
values 5* are shown versus J 3 . The solid (dashed) lines de- 
scribe the orthogonal (unitary) symmetry. Inset: the stan- 
dard deviation of the peak spacing cr(Aa) versus J s for the 
even, odd and combined ("total") transitions. J s and cr(A2) 
are measured in units of A. 

Koopmans' limit is exact when the fluctuations of the 
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matrix elements are ignored, and our calculations of A2 
reduce to those in Ref. [O, i.e., they can be derived di- 
rectly from an Hamiltonian that consists of a random 
one-body part plus an average interaction (0). In the 
limit g — > 00 this Hamiltonian is just the universal Hamil- 
tonian with J s = A/2 [ p^,p^[ . A better estimate of J s 
can be obtained in RPA flllf; it increases monotonically 
with r s but remains below A/2. 

In the simple limit (R), the spin and spacing distribu- 
tions are determined by a single parameter J s /A [pj. 
This limit is demonstrated in Fig. |l| where we show the 
probabilities of various spin values versus J s . The inset 
shows the standard deviation of A2 versus J s for the even 
and odd transitions as well as the combined one ( "total" ) . 

Next we discuss the fluctuations of the interaction ma- 
trix elements |f2(i| , p4| , which are approximately Gaussian 
variables. We discuss separately the bulk screened inter- 
action and surface-charge potential. Using the diagram- 
matic approach for the two-body screened interaction, we 
have, for r s <C 1 and to leading order in 1/g 



GOE : a(v a p) 
GUE : a(v aP ) 
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pj — v ^02; <J(v c 
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= 2V2a 2 
V2a 2 ■ (8) 



Different matrix elements (including the direct v a p and 
exchange v^g) are uncorrelated. We note that the coef- 
ficients of ct 1 and 02 in Eqs. are different from those 
obtained for a zero-range interaction |25f| . The parameter 
(T2 is given by 



02 



X 



dr 2 U 2 {r 1 ,r 2 ) 



1/2 



c 2 A/ 5 , (9) 



where c 2 is a geometry-dependent coefficient. For a disk 
of radius R and boundary conditions of vanishing nor- 
mal derivative, we find g!§ c 2 = 2E^ m ^] 1/2 « 0.67 
where x^ m are the zeros of J[ (x) (Ji is the Bessel function 
of order I) and gA — 2ir T1D/R 2 (D being the diffusion 
constant). Since only a few matrix elements contribute 
to the peak spacing, the contribution of the two-body 
screened interaction is parametrically of the order A/g, 
unlike the A/y/g dependence found in Ref. |l2] ]. 

The surface-charge contribution to an interaction ma- 
trix element is v aj 3 = V a +Vp, where V a = J \i/j a (r)\ 2 V(r) 
is a diagonal matrix element of the surface-charge poten- 
tial. We have 

a(V a ) = (2/p) 1 / 2 a 1 ; T£T£ - T£ T£ r» , (10) 

where 

1/2 



A' 2 J dri j dr 2 V(r 1 )n(r 1 ,r 2 )V(r 2 ) 
C1A/V5 ■ 



(11) 



For an isolated two-dimensional (2D) circular disk of 
radius R, the surface-charge potential can be approx- 



k = 2ire 2 is/e is the inverse screening length in 2D and 
e is the dielectric constant EG]. We then find c\ = 
2" 1/2 E m# oSin 2 x 0im /(4 m jJ(i 0im ))] 1 /2 „ 0.087 Q. 

The above results can be generalized to a ballistic 
dot, using the ballistic supersymmetric a model obtained 
when a weak disorder with finite correlation length is 
added |2(|. In particular, if the Lyapunov length of this 
smooth disorder is smaller than the dot's size, relations 
similar to Eqs. (||) and ( |l0| ) can be derived but with the 
ballistic propagator replacing the diffusive propaga- 
tor n in Eqs. ([)]) and ([ll]). For a circular dot we define 
the ballistic Thouless conductance from the inverse time 
it takes to cross the diameter 2R of the dot, leading to 
g = nkpR/2 = 7r(n/2) 1 / 2 (n is the number of electrons 
in the dot). In the ballistic case 

(72 = c 2 A[\n(c' 2 g)]^ 2 /g , (12) 

where c 2 = y3/2 is a geometry-independent constant. 
The direct propagation between ri and r 2 contributes 
to Hb('T'i,i'2) a term l/(irkp\ri — r 2 \) which at shorter 
distances should be replaced by its quantum counter- 
part Jo(^H r i — r 2|) @- The corresponding contribu- 
tion needs to be renormalized such that its integral over 
T\ (or r 2 ) vanishes J28|. The remaining part of in- 
volves single or multiple scattering from the boundaries 
and calculating it requires knowledge of the semiclassical 
dynamics. It can be calculated analytically for a circu- 
lar billiard with diffusive boundary scattering [Q . Using 
this model, we estimate c' 2 = 0.81 (lj]]. A similar estimate 
of @) gives ci = 0.123. 

We have studied the effects of fluctuations of the inter- 
action matrix elements on the peak-spacing distribution. 
In Fig. ^| we show the standard deviation c(A2) versus 
(7i ( J s = 0.3A and a 2 = 0.05A) for the orthogonal (solid 
lines) and unitary (dashed lines) symmetries. cr(A2) in- 
creases with (7 1 and it does so faster in the odd case. 
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FIG. 2. The standard deviation <r(A2) versus the standard 
deviation o\ of the surface-charge potential for J s = 0.3A 
and 02 = 0.05A. The solid and dashed lines describe the 
orthogonal and unitary symmetries, respectively. 

Fig. |3| shows typical peak-spacing distributions for 
both the orthogonal (left) and unitary (right) symmetries 



3 



for cr 2 = 0.025A and cr x = 0, 0.03A and 0.06A. Signa- 
tures of the bimodality can still be observed for ax — 
but they disappear at a% = 0.03A. Nevertheless, the 
distributions remain asymmetric (more so in the unitary 
case). 

An additional contribution to the fluctuations of A2 
arises from the variation of the gate voltage between 
peaks. In general, the change of the gate voltage be- 
tween two peaks leads to a spatially non-uniform change 
U(r) = — V(r) + V(r) in the confining potential, where V 
originates in the mutual dot-gate capacitance |l5[ ]. This 
leads to scrambling of the HF levels between peaks. For 
example, let us consider the odd transition. As the gate 
voltage changes between V™ +1 and Vg l+2 , the reference 

HF levels e£° change by Set ] w U a = J dr\iP a (r)\ 2 U{r). 
We therefore substitute — > + U a in the calcu- 
lation of /i n +2 and in Eq. (^) {fi n +i is unchanged and 
calculated from Eq. (||)). This level scrambling can also 
lead to spin rearrangement in the dot. The ground-state 
spin of the n + 1-electron dot at gate voltage Vg +2 (just 
below the transition) is found from (^) (and its gener- 
alization to higher values of S) after the substitution 



» 



+ U a . 
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FIG. 3. Peak-spacing distributions P(A2) for the or- 
thogonal (left) and unitary (right) symmetries and for 
01 = 0,0.03A and 0.06A. In all cases J s = 0.3A and 
(72 = 0.025A. The bimodality is lost already for a\ = 0.03A 
but the distributions remain asymmetric. 

The fluctuation properties of U a are similar to those 
of V a in Eqs. (|Io|), except that V(r) is replaced by U(r) 
in Eq. ©. 

We now compare our theory with the experimental re- 
sults of Ref. H at the lowest measured temperature of 
T = 0.22A. At this temperature it is necessary to include 



the effect of excited states and in particular the contri- 
bution from both lowest S = and S — 1 states [ p^| . 
We model the gate- voltage scrambling by a potential V 
whose matrix elements are uncorrelated from the matrix 
elements of V and have the same variance. It is difficult 
to estimate o\ and 01 for the ballistic dot used in the ex- 
periment. The simple estimates based on a billiard with 
diffusive surface scattering (see the second paragraph af- 
ter Eq. (pd])) for n w 340 electrons (i.e., g ps 41) give 
(7i = 0.02A and (72 = 0.04A. Fig. [I] compares the exper- 
imental distribution of Ref. [|| in the presence of a mag- 
netic field (solid histograms) with the corresponding the- 
oretical distribution (dashed histograms) that includes an 
experimental noise of 0.1 A. The theoretical distribution 
describes rather well the asymmetry of the experimental 
distribution and its width 0.27 A is slightly below the ex- 
perimental width of (0.29 ± 0.03)A. In chaotic billiards 
our estimates of ui and 172 can be enhanced by up to a 
factor of 2 and lead to better agreement with the data. 




-0.5 0.5 

A 2 -<A 2 > 

FIG. 4. Calculated peak-spacing distribution at T = 0.22A 
for the unitary symmetry (dashed histograms) and 
J s = 0.28A, (7i = 0.02A and a 2 = 0.04A is compared with 
the experimental distribution of Ref. || (solid histograms) 

In conclusion, we have developed a HF-Koopmans ap- 
proach to study spin and interaction effects in diffusive 
or chaotic quantum dots. In particular we have studied 
the dependence of the peak-spacing distribution on the 
fluctuations of the interaction matrix elements to leading 
order in the inverse Thouless conductance. We find good 
agreement with the lowest temperature data of Ref. || . 
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